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Abstract - We report an unusual giant linear magnetostrictive effect in the ferrimagnet 
Gd 2 /3Ca 1/ / 3 Mn03 (T c «80 K). Remarkably, the magnetostriction, negative at high temperature 
(T « T c ), becomes positive below 15 K when the magnetization of the Gd sublattice overcomes 
the magnetization of the Mn sublattice. A rather simple model where the magnetic energy com- 
petes against the elastic energy gives a good account of the observed results and confirms that 
Gd plays a crucial role in this unusual observation. Unlike previous works in manganites where 
only striction associated with 3d Mn orbitals is considered, our results show that the lanthanide 
4/ orbitals related striction can be very important too and it cannot be disregarded. 



Introduction. — Manganites are perovskites mostly 
known for their spectacular colossal magnetoresistance 
(CMR): the electrical resistivity can change several orders 
of magnitude under a moderate applied magnetic field B 
[1] . They also show another impressive property called gi- 
ant linear magnetostriction (MS): sample dimensions are 
strongly affected by a magnetic field, either external or 
molecular [2]. The effect is comparable in magnitude 
(AL/LMCT 3 at several Tcsla) to the highest MS values 
ever reported. Both CMR and MS arc particularly large 
around the Mn-ions ferromagnetic ordering temperature 
[3]. Associated with this order, and depending on the dop- 
ing level, manganites can display a metal-insulator (MI) 
transition, too. In this way, manganites offer a unique 
testing ground to study the interplay between electronic, 
spin and lattice degrees of freedom. 

As expected, structural distortion of the plain per- 
ovskite structure strongly affects the magnetic and 
electronic properties of manganites. This is usually 
parametrized by the so called tolerance factor t, which 
quantizes the mismatch between the size of the different 
ions in the formula. This mismatch primarily influences 
the exchange interaction between Mn ions altering both 
the length and the angle of the Mn-O-Mn bond. An 'uni- 
versal' temperature T versus t = {d R /A-o) I V2(d/un-o) 
phase diagram has long been reported [4] for the hole 
doped R 2 /3A 1 / 3 Mn03 manganites (R is a lanthanide and 



A is an alkaline-earth element). Slightly distorted struc- 
tures (t ~1) show the insulating-paramagnet (PMI) to 
metallic-fcrromagnet (FMM) transition. However, the 
metallic state disappears at higher distortions (£<0.91) 
even though a transition to a insulating-ferromagnet 
(FMI) is observed. 

An estimated value of tw0.89 places Gd 2 /3Ca 1 / 3 Mn0 3 
well inside the insulating regime. Indeed, no MI transition 
is observed down to 5 K with the resistivity p showing a 
characteristic semiconducting behavior in the whole tem- 
perature range 016]. Nevertheless, magnetic properties 
are quite more interesting. Mn magnetic moments start 
ordering ferromagnetically around T c ~80 K. The Gd mo- 
ments react to the internal field created by the Mn fer- 
romagnetic sublattice gradually aligning in the opposite 
direction. The two sublattices compete each other giving 
rise first to a maximum in the magnetization around 50 K 
and finally to a full compensation at T com p~15 K where 
the magnetization vanishes. At lower temperature, the Gd 
magnetic moment overcomes the Mn moment. The overall 
temperature dependence of the magnetization corresponds 
then to a ferrimagnet created by the two opposite Mn and 
Gd sublattices [HErE]. 

In this work we study the magnctostructural proper- 
ties of Gd2/3Cai/3Mn03. The rather complex magnetic 
structure clearly couples to the atomic lattice giving rise 
to a giant linear magnetostrictive effect [10] . Remarkably, 
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the negative field dependence of the MS observed at high 
temperature changes its sign and becomes positive when 
T < T comp . We use a 4-site mean field approximation to 
model the experimental data. It demonstrates that the 
competition between Gd-Gd and Gd-Mn spin correlations 
is responsible of the sign change in the MS. This finding 
shows that the usually underestimated MS associated with 
the lanthanide Af orbitals in manganitcs can be compara- 
ble to the usual giant striction given by the re-orientation 
of the Mn 3d orbitals. 

Experimental details. — Pure single crystalline 
samples of Gd 2 /3Ca 1 / 3 Mn03 were grown by the floating 
zone technique. Crystal quality and composition have 
been checked through XRD and EDS scans. A capaci- 
tive technique was used in the dilation experiments. The 
high resolution (<1 A) dilatometer [TT] is placed in a evac- 
uated environment with a low pressure (P <10 _1 torr) of 
exchange He 4 gas. Magnetic field is applied along the [020] 
direction of the orthorhombic Prima crystalline structure 
(a= 5.39A, 6= 5.56Aand c= 7.5A) in all the experiments. 
Dilation experiments are always performed in a longitudi- 
nal configuration with B \\ L | [020]. Several samples of 
different sizes have been measured with a perfect agree- 
ment between them. Sample length L is typically about 
200 fxm. 

Results and discussion. — Representative isother- 
mal linear magnetostriction results after a zero field cool- 
ing procedure are shown in Fig. [T] (solid lines). The effect 
is giant with no evidence of saturation up to B = 12 T (the 
highest applied field), reaching a maximum value around 
T c -80 K (AL/L «10~ 3 ). Two very distinctive regimes 
are found: 

(i) above T comp ~15 K the field dependence of L is neg- 
ative and monotonic. Hysteresis and relaxation effects arc 
important, mainly in the range 40 K< T <T C , 

(ii) below T cornp magnetostriction becomes positive at 
low fields B < 7 T (the initial negative slope at B < 1 T 
is associated to magnetic domains and is absent if the ex- 
periment is performed after a field cooling procedure). At 
higher fields it turns negative again resulting in an overall 
non- monotonic field dependence of the MS. On the other 
hand, around T comp (where the magnetization almost van- 
ishes), the magnetostriction is negligible below B^ 4 T. 

Tcomp marks the onset of the Gd magnetic ordering 
which dominates the low temperature regime while Mn 
moments prevail in the high temperature regime. In this 
sense, this unusual magnetostriction strongly points to- 
ward the interplay of the different magnetic interactions: 
Mn-Mn, Mn-Gd and Gd-Gd. We use a simple model to 
verify this hypothesis where the three different interac- 
tions are introduced in the Hamiltonian via Hcisenbcrg- 
like terms. 

We consider an homogeneous network of Gd ions with 
one Gd ion for each Mn, ignoring the random nature of 
their localizations. As there are 2/3 Gd ions for each Mn, 
we rescale the Gd effective magnetic moment to J= 2/3 x 
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Fig. 1: (color online) Experimental (solid) and calculated 
(dashed) magnetostriction. Upper (lower) panel shows results 
in the low (high) temperature range. Curves are vertically 
shifted. 

7/2 = 7/3 ~ 5/2. We choose the 5/2 value for the Gd spin 
to retain the quantum nature of the spin without rescaling 
of the g- factor (gad = 2). Manganese ions appears in a 
mixture of 1/3 of S=3/2 and 2/3 of S=2. As in both 
cases the orbital magnetic moment is quenched we take 
<?Mn = 2. To keep the experimental zero temperature net 
magnetic moment of 1 fj.B (perfect ferrimagnetic ordering 
given by the two sublattices [5]) we take S=2 for the Mn 
effective magnetic moment. 

Based on these considerations we write a Hamilto- 
nian with a ferromagnetic (coupling K Mn-Mn) network 
of manganese spins S (S=2) antiferromagnetically cou- 
pled (KMn-Gd) to a network of ferromagnetic (Kod-Gd) 
gadolinium spins J (J— 5/2). The magnetic interaction is 
given by the Hamiltonian: 

H m = K Mn _ Mn ^ ■ Sj + K Mn _ Gd 

(id) » 

+K G d-Gd Ji ■ Jj + 9HbB ■ ^(Si + Ji) (1) 

(id) 1 

The smaller values of the effective spins allow us also to 
use a 4-sites (2 Mn and 2 Gd) cluster in the Constant Cou- 
pling approximation (see Appendix) [T2"Hl"4"] . We consider 
six neighbours (2 = 6). In the Constant Coupling approx- 
imation the interactions are isotropic, meaning that they 
represent averaged interactions. 
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Following early works [T4TU6j . we consider that the ex- 
change parameters are strain dependent. If the lattice is 
under some small distortion 8l all the coupling parameters 
change accordingly: 



(2) 



From fittings to magnetization experiments, we obtain: 
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K and K Gd _ Gd = K. 



Mn-Mn coupling K^ /In _ Mn is ferromagnetic with T c close 
to the experimental value of ~ 80 K; Kj^ Gd is antiferro- 
magnctic and gives T comp ~ 15 K together with an effec- 
tive null coupling between Gd ions (also expected due to 
the dipolar origin of those interactions) . Both, our exper- 
imental and calculated magnetization results are similar 
to those reported by Snyder et al. [5]. 

The presence of a distortion also increases the elastic 
energy 

E e = 1/2C6 2 L (3) 



For a state \G), the total energy E m + E e is minimized 



for 



6l(B) = -^(G\a^Si-Sj+fi'^Si-Ji+i'Y^Ji-Jj\G)B 

As we are interested in the length distortion respect to 
the B = case, we compute 
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(G\Y,S l -S ] \G) B ) 

(id) 

(G^Si-JiMs) (4) 
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where 



h and x 



{a, (3, 7}. ()b denotes the thermal 



y x ~ c 
expectation value. 

The correlations are computed using the 4-site Constant 
Coupling approximation and correspondingly now i,j = 
1,2. Figure [2] shows the computed spin correlators as a 
function of magnetic field (A(O) = (O) b — (O) b=o, where 
O = Si ■ S2, Si ■ 3\ or J\ ■ J2) in the different temperature 
ranges: below, around and above T comp . 

There are several issues to emphasize: 

i) the field dependence of the Mn-Mn correlations 
(A(Si • S2)) is monotonic and, except around T c , it is 
also very small (Fig. HJa)), reflecting the fact that Mn 
sublattice is almost fully polarized at low temperature. 

ii) Mn-Gd correlations (A (Si ■ Ji)) show the largest ef- 
fect below T c , as seen in Fig. [U[b). This is reasonable 
since the applied field tends to align both sublattices in the 
same direction gradually destroying the otherwise almost 
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Fig. 2: (color online) Computed correlations change respect 
to the zero field situation for 8 K (below T comp ), 14 K (close 
to T comp ) and 45 K (above T comp ). (a) Mn-Mn correlations, 
(b)Mn-Gd correlations and (c) Gd-Gd correlations. 



perfect low temperature ferrimagnet. Its contribution to 
the low field total magnetostriction, however, is smaller 
than that due to Gd-Gd correlations by a factor j approx- 
imately. At higher fields, on the other hand, A(Si ■ J2) 
becomes the larger contribution. 

iii) Gd-Gd correlations (A(Ji ■ J2)) show the most rele- 
vant and distinct behavior. While A(Ji ■ J2) is negative for 
T > T comp . it vanishes around T comp and becomes positive 
when T < T comp . The physical origin of this sign change 
is related to the local field acting on Gd moments. This 
local field is made up of the molecular field created by the 
Mn moment Bmu plus the external field B. Above T comp , 
Bun points opposite to B so an increasing B (B < Bmu) 
results in a decreasing local field. This decreasing local 
field reduces the Gd moment and consequently the Gd-Gd 
correlations. Below T comp , on the other hand, Bmu points 
in the same direction of B, so an increasing B (B < Bmu) 
results in a increasing local field. This increasing local field 
raises Gd-Gd correlations. This sign change in the Gd-Gd 
correlations is indeed responsible for the sign change in 
the magnetostriction observed at low temperature. 

Computed MS curves are also shown in Figure [TJ (A Q = 
12 x 1(T 4 , Ay3 = 0.8 x 10~ 4 and A 7 = -1.6 x 10~ 4 ). A Q 
is chosen as to get a good agreement at high temperature 
(T > T c ) where the only non-negligible correlator is A(Si ■ 
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S 2 ) while A 7 is chosen as to get a positive striction at low 
temperature and field (T < T comp , B < 5 T) where only 
A(Ji • J 2 ) is non-negligible. is then selected to get the 
best agreement in the whole temperature and field range. 

The model gives a good account for the non-monotonic 
AL/L below T comp (see curve at 8 K, Fig. QJa)). It is 
a consequence of two opposite contributions: a negative 
Gd-Gd magnetostructural coupling (A 7 < 0) and a pos- 
itive Mn-Gd coupling (A^ > 0). As stated previously, 
Mn-Mn correlations are almost saturated and they do not 
contribute to the MS in this low temperature range. The 
coupling parameters (Ks) used (fixed by the fit of the mag- 
netic properties) allows us to reproduce only qualitatively 
the field value of the MS maximum. The model also ac- 
counts for the extinction of this maximum at T comp ^15 
K, where the striction becomes very small. 

At higher T {T comp < T < T c ) the MS gets negative. 
As before, Mn-Mn correlations almost do not change but 
now both Gd-Mn and Gd-Gd correlation effects point in 
the same direction. In the intermediate and high tempera- 
ture range (T T comp ), where the Mn-Mn contribution is 
the more relevant one, the agreement is fairly good, except 
around T c . Not only the magnitude of the magnetostric- 
tion is well accounted, also the curvature of the isotherms 
is very well reproduced (see Fig. QJb)). It is interesting 
to stress that even though magnetization is an increasing 
function of field in the whole temperature range, the mag- 
netic correlations and so magnetostriction shows two very 
distinctive regimes: a monotonic magnetostriction at high 
temperature that becomes non-monotonic below T comp . 

In this isotropic model, C can be estimated as vBt, 
where v is the volume of the perovskite unit cell and Bt 
is the Bulk modulus. Mn-Mn parameter A Q is positive, 
and so is a. Since the interaction is FM {Kun-Mn < 
0), that implies that \KMn-Mn\ increases as the lattice 
gets smaller. This is the expected behavior for exchange- 
like interactions. There are no available pressure effects 
on Gd 2 /3Ca 1 / 3 Mn03 to compare with. Nevertheless, a 
rough estimate can be done by replacing Gd by another 
lanthanide, i.e. by chemical pressure. This substitution 
(keeping the composition at R 2 /3Cai/ 3 Mn03) does not 
modify the Mn valence and so, its magnetic moment. So, 
a change in T c can in principle be associated with an inter- 
ion distance d change. In this isotropic approximation, 
d = v 1 / 3 . Taking B T = 150 Gpa [TZ], and v = 55.36A 3 , p] 
we get a = 722 K. For Dy 2/ 3Ca 1/3 Mn03, v Dy « 55.10A 3 
[T5] . This results in a change of the exchange parameter 

d-Gd 



K 



Mn—M 



n given by AK 



-2 K. This very 
rough estimate of a 20 percent increase in the exchange 
parameter [KMn-Mn = -9 K) is of the same order of mag- 
nitude that the T c increase observed in Dy 2 /3Cai/3Mn03 

Gd-Gd parameter A 7 is negative and it may be related 
to the dipolar (anisotropic) origin of these interactions: 
for zero distortion the net (average) interaction is zero but 
when the lattice shrinks AF interactions prevails. Notably, 



the magnetostriction associated with Gd (low field positive 
striction below T comp ) is of the same order of magnitude 
than the observed striction in metallic gadolinium |14j . 
On the other hand, Mn-Gd parameter Ap is positive and 
it is much more difficult to understand since Mn-Gd in- 
teractions are antifcrromagnctic: this effective interaction 
diminishes as the lattice shrinks. This counter-intuitive 
value of A^ could be related with the rotation of the oxy- 
gen tetrahedra that sourrounds Mn ions. 

Conclusions. — An unusual non-monotonic giant 
magnetostriction is observed in single crystals of the fer- 
rimagnet Gd^CajysMnOs at low temperature (T < 
Tcomp ~15 K) arising from the interplay between the Mn 
and Gd magnetic sublattices. A simple mean field ap- 
proximation where different magnetic interactions (Gd- 
Gd, Mn-Mn and Mn-Gd) compete among them and with 
the elastic energy gives a good account of the observed 
results. Particularly, the change in the sign of the magne- 
tostrictive effect at low temperature is driven by the com- 
petition between Mn-Gd and Gd-Gd magnetic correlations 
and it does not involve the Mn-Mn correlations. Unlike 
previous works |20j in manganites where only striction as- 
sociated with d orbitals is considered, our results show 
that / orbitals related striction can be as important. 

Appendix. 

Constant Coupling Approximation . The constant 
coupling (CC) approximation [T2Trr4] is an improvement 
over a classical mean field approximation. It allows to 
consider correlations and to obtain a critical temperature 
closer to the exact result. A classical mean field approxi- 
mation replaces all the interactions of a site by an effective 
magnetic field made up of the external field and its neigh- 
bours magnetization. This neighbours magnetization is 
assumed to be the same than that of the site. In the con- 
stant coupling approximation two systems must give iden- 
tical results for the magnetization. If the original problem 
is in a network with z neighbours per site, one system is 
made with a single site and z "effective" neighbours, while 
the other is made up with a cluster of 2 sites and z — 1 
"effective" neighbours for each site. In a classical mean 
field we have to search for an effective field proportional 
to the neighbours magnetization. The CC approximation 
consists in searching for an effective field in both systems 
such that the same magnetization is obtained in the single 
site and the cluster. 

To illustrate the procedure we use a simple spin 1/2 
fcrromagnet. For a single site the Hamiltonian is written 
as 

H ss = g^B(zB eff + B) ■ S 

where i? e // is an effective field. In a classical mean field 
approximation i? e // = KMn-Mn{S) and a self-consistent 
(S) is looked for. For a two sites cluster the Hamiltonian 



H 2s = KSx ■ S 2 + gnu (z-l)B eff + B ■[S 1 + S 2 
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Fig. 3: (color online) Magnetization as a function of the ef- 
fective field for a single site and a two site cluster (S=l/2, 
K — 30K, z — 6). Upper panel: 27 K. Lower panel: 20 K. 
Insets show the magnetization difference. The classical mean 
field transition temperature is Tc,mf = 2/3zS(S + 1)K = 90 
K. 



In the upper panel of Fig. [3]we show the magnetization 
at high temperature (27 K, just above the transition) of 
both a single site and a cluster of two sites as a function of 
B e ff with an applied external field B = 0.1 T. These mag- 
netizations are zero for (z—l)B e ff — —B (or zB e ff = —B 
for a single site) and show the expected paramagnetic-like 
behaviour for small systems. At temperatures above the 
transition temperature both magnetizations agree for a 
single field which is the searched effective field (around 
7 T for this temperature; the corresponding magnetiza- 
tion is less than 0.1/is). As the temperatures lowers, this 
solution moves toward higher fields with a corresponding 
larger magnetization. Below the transition temperature, 
two new solutions (higher in energy) appears, just as in a 
regular mean field approximation (lower panel of Fig. [3]) . 

For the ferrimagnet Gd 2 /3Ca 1 / 3 Mn03 we take as the 
"single" site an unit made up of an effective manganese 
(S=2, g = 2) and an effective Gd (J=5/2, g = 2) ion 

H ss = gLiB(zB e ff,Mn + B)-S 
+ gfi B (zB ef f tGd + B) ■ J 

+ Kad-MnS ■ J 



The cluster is made with two Mn and two Gd ions and 
we take z = 6. Each Mn (Gd) ion interacts with the other 
and with an external field made by the z — 1 remaining 
neighbours. In each site, there is a Gd-Mn interaction. 



The Hamiltonian is 

H-2s = KjUn-MnSl • S2 + Kcd-Gd'h ' 'h 
+ K Gd-Mn{S\ ■ J\ + S2 ■ J 2) 

+ gfi B {z - l)B e ff t Mn + B ■ [Si + S 2 ] 
+ guB (z - l)B eff . Gd + B • [ Ji + J 2 ] 

When the effective field is found, correlations between 
the different ions can be computed in the 2-site (four ions) 
cluster. 
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